Noise in superconductor- quantum dot-normal metal structures in the Kondo regime 
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We consider a N-dot-S junction in the Kondo regime in the Hmit where the superconducting gap is 
much smaller than the Kondo temperature. A generalization of the floating of the Kondo resonance 
is proposed and many body corrections to the average subgap current are calculated. The zero 
frequency noise is computed and the Fano factor sticks to the value 10/3 for all voltages below the 
gap. Implications for finite frequency noise are briefly discussed. 
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Possible realizations of quantum dots have revived the 
interest in Kondo physicsii^. For normal electrodes, 
the dot spin is screened and the system behaves in 
first approximation as a resonant level at the Fermi en- 
ergy. However, this picture is not sufficient; the sys- 
tem is a Fermi liquid with interactions, morevover, the 
ratio of elastic and inelastic backscattering obeys some 
universality'^, which can be obtained from the concept 
of floating of the Kondo resonance. Current noise is a 
useful tool to obtain informations on interactions, which 
are not present in the average current (/) . Recently, the 
zero frequency noise S and the Fano factor F — S/2e{I) 
in the SU(2) case^ and the SU(4) case-'^ were calculated 
and confirmed experimentally^ for SU(2). 

What happens to this effect in the superconducting 
case? This poses the problem of interplay between 
Kondo physics and superconductivity, already present 
in heavy fermion compounds^ and underdoped high-Tc 
materials^. Devices with a dot between two supercon- 
ducting electrodes have been extensively studied both 
theoreticallySiiS'iiii^'i^ii^'i^ and experimentallyi^iii with 
emphasis on the Josephson current, and on how increas- 
ing the gap A destroys the Kondo effect. Here, we study 
the noise of the subgap Andreev current in a normal 
metal-dot-superconductor structure for A ^ Tk- While 
in this regime, there is still no destruction of the Kondo 
resonance by the presence of the gap Afii^l, the inter- 
play between one and many-particle scattering, Andreev 
and normal reflexion, is far from trivial. The central re- 
sult of this paper is a generalization of the floating of 
the Kondo resonance in the case where one electrode is 
superconducting. The most noticeable consequence is a 
constant Fano factor, equal to 10/3 for all voltages bias 
below the gap. This could be tested experimentally on 
carbon nanotubes. 

Model: A quantum dot with effectively one level of 
energy eo is placed between two electrodes, the left one 
being normal and the right one being a usual BCS super- 
conductor, see Fig. 1. The on site repulsion U on the dot 
is supposed to be the largest energy of the problem. Elec- 
trons can hop from the lead to the dot with amplitude t, 
implying a broadening of the level F = 2Trp{€F)\T\^ , with 
piep) the density of states at the Fermi energy in the 



normal metal. We abide by the particle-hole symmetric 
case, for which eo = — U/2. The Hamiltonian reads 

H = ^ (ffc - A*p)4,<T,pCfc,<T,p + TCkM,pd]j + h.c. 
fc 

where A is the superconducting order parameter; p = L 
for the normal lead on the left and p = i? for the su- 
perconducting one; the operators cj, ^ ^ and dl create an 
electron on the p lead and on the dot, respectively and 
Tier — dl^d^- The chemical potentials are fi^ = and 
PR = 0, with V the voltage bias. For A ^ Tk, and even 
when both electrodes are superconducting, it was shown 
by various methods^, both analytical and numerical^^'— 
that the Kondo resonance subsists despite the gap. In 
the non-Kondo case, the average current has been calcu- 
lated in Ref. [l^. Also, in the low temperature regime, 
in the Kondo case, the average current and the noise 
were estimated by using slave-boson methods but this 
method is an effective one-body approximation to the full 
Kondo Hamiltonian and thus, is not sufficient to capture 
the complexity of the Kondo Hamiltonian, already when 
both electrodes are normal''. 

Method and one-body setting: Here, we want to gener- 
alize the calculation of the Kondo noise when two elec- 
trodes are normal to the case where one electrode is su- 
perconducting. In order to include the many-body ef- 
fects, we consider a slightly different model for the Kondo 
dot, which has been used for modelling an imperfect NS 
junctiorti^, see Fig. 2. The dot is a scatterer placed 
in front of the superconductor and produces dephasing, 
whereas in reality, the dot acquires off-diagonal order due 
to hopping on and from the superconductor. An incident 
electron of energy e is scattered by the dot and then An- 
dreev reflected as a hole by the NS interface, supposed to 
be perfect, with a modulus one and a phase —i atan(e/A), 
and then, this hole gets scattered by the dot, (see Fig. 
2). The dephasing of the electron by the scatterer is 

= Si°^ + Si^^ + where si°^ = -7r/2, si^^ is first 
order in e/Tx and the dots represent higher orders in 
e/Tff. The energy transmission of electrons through the 
scatterer is T(e) = sin^((5e). The same thing happens 
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FIG. 1: Quantum dot sandwiched between a normal electrode 
on the left with chemical potential /^l an a superconducting 
one on the right with gap A and chemical potential fiR — 0. 
Only one level of negative energy eo intervenes. Hopping am- 
plitudes on and from the dot to both electrodes are assumed 
to be equal. The effective energy level taking into account the 
Hartree correction is eo = eo + U/2, in heavy dashed line. 



for the holes with dephasing Sh- Using the Bogolubov-de 
Gennes (BdG) equations enables to express the s— matrix 
in Nambu space for the whole structure, which is 4 x 4 
but, below the gap, reduces to a 2 x 2 reflexion matrix, 
because no transmission of one-particle excitation occurs 
in the superconductor. The energy dependent normal 
reflection coefficient and the Andreev amplitude re- 
flection coefficient r ^ can be extracted in terms of Se and 
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FIG. 2: Simplified model for the A'' — dot — S system. A scat- 
terer, having no off-diagonal Green's function, is placed in 
front of the superconductor accounts for the imperfect trans- 
mission properties of the dot. The interface with the su- 
perconducting electrode is taken to be perfect so that only 
Andreev reflexion occurs at the interface and only normal 
reflexion or transmission occurs across the scatterer. Prox- 
imity effects are neglected. An incident electron is partially 
transmitted through the scatterer and then totally Andreev 
reflected as a hole. This hole is in turn partially transmitted 
as a hole but also partially reflected by the scatterer, etc... All 
processes need to be resummed to infinite order, the ampli- 
tudes rjv and va are the result of this (coherent) resummation. 
However, this is equivalent to solving the BdG equations. 

Sh- We want to make contact with the original model of 
Eq (I). If Tr- were infinite, then, the one-body picture of 
a resonant model centered at the Fermi level, with width 
Tk would be valid. We are close to the unitary limit 
and we need the dot Green's function in this limit. In 



the case where both electrodes are normal, the retarded 
Green's function is of the form z"^ (e — eo — 2iTsgn{€)) 

with z^l- (aS'^'=*(e)/9e)|,=ep and f = z'^T ~ Tk/2; 
S'''^*(e) being the retarded self-energy of the dot and 
eo = Co + U/2. This has been used justified theoreti- 
cally by Ref. d When both electrodes are supercon- 
ducting, for A <C Tfc, a similar procedure exists and has 
been used in Ref. [ij to study the dynamics of Andreev 
states. We thus use these Green's functions, adapted to 
our case, and following the same steps as in Ref. [lO|, the 
conductance gA reads 
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On the other hand, the model of Ref. \2M leads to a con- 
ductance 
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Thus, for gA and gB to have the same expression, we ad- 
just the dephasing of the scattering center to be — 
(e/A) (e/r) = (e/Tft-) ai(e), where ai is a function of e. 
The same occurs for 5^^^^ . This form does not correspond 
to a resonant level at e = and width r/2. In the original 
model, the dot acquires some off-diagonal matrix element 
in Nambu space. Very qualitatively, this favors direct An- 
dreev reflexion of a wave incoming on the dot. Since, for 
the whole structure, |rArp + |rAp = 1, there will be less 
normal scattering. The amount of normal scattering in- 
creases as {e — eF)"^ instead of [e — ep) for the normal case. 
Let us denote by G^q^ i i the upper Nambu component 
of the retarded dot Green's function, calculated with- 
out the many-body corrections, using the renormalized 
parameter F and putting aside the usual multiplicative 
renormalization factor z. For energies higher than the 
gap, the form of the normal part of the spectral density 
of the dot pi,i(e) = — tt"^ ^"^G^oti 1 '^'^'^y close to the 
usual shape for the normal case. Surprisingly, this per- 
sists down to energies barely larger than A. Of course, 
this is not the exact pi,i(e) but this feature seems to 
be shared by more elaborate numerical solutions for the 
S-dot-S case, such as numerical renormalization group 
(NRG)^^ or functional renormalization group (fRG) ^"^'^^ . 

Many-body calculation: Now, we want to put the 
many-body corrections since Tk is not infinite. The elec- 
trons build the Kondo resonance from hopping to the 
normal lead and also from hopping to the superconduc- 
tor. Unlike the normal case, the argument of doping 
the normal side to establish the floating of the Kondo 
resonance^^ does not work directly, because doping the 
metal and changing the chemical potential for the elec- 
trons will also change it for the holes. However, for 
A ^ e ^ Tk, 5^e^ is linear in e, (ai becomes energy 
independent) and the argument can be used. Thus we 
have to cancel the linear contribution in 1/Tk of Se by 
a many-body interaction between electrons (and holes). 
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as in usual Fermi-liquid theory. The fixed point Hamil- 
tonian is 



Hpp — Hq + Hint , 



(4) 



with Hq being the one-body Hamiltonian, involving scat- 
tering states (first order in l/Tj^) and Hint is the inter- 
action between quasiparticles. The simplest form main- 
taining particle-hole symmetry is 
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with bk 

J 



,<T,i?)/\/2 and Bk,a = {cL^.w , 



-k.Tj,R 



)/\/2 with 'a — — (T and /3i(e) is a function of en- 
ergy. Here, a right-mover state with energy e will have 
some normal reflected part. No elastic scattering Hamil- 
tonian is necessary^. 

Because of Hartree-Fock corrections to order one in (3i , 
the dephasing s'^p will be changed to (5e^^ — {l3i/TK)5n-ff, 
so that (3i = ai to cancel the 1/Tk overall contribution. 

The main assumption of this paper is that this can be 
continued all the way down to e = 0. Then, /?i(e) will 
have to cancel ai{e). 

Now, we calculate the average current in the Keldysh 
formalism, perturbatively to order following the 

method of Rcf. 5. The right mover scattering states are 
thus of the form, in Nambu space, far on the left of the 
scattering center 



^n{ek,x) = 



(e +rive 



— ika 



(6) 



The left-mover scattering state consists in sending a hole 
from the left, which is partially normally reflected and 
also Andreev reflected as a left propagating electron. At 
the one-body level, the current operator has two compo- 
nents and is then expressed in terms of the BdG wave 
function, and falls into three parts, Ir, involving only 
right-moving scattering states, II with solely left-movers 
and loD which involves both right and left- movers. Each 
of these three terms gives rise to two contributions. One 
contribution, denoted by subscript 1 comes from the elec- 
tron part of the two component Nambu wave function 
and the other one from the hole part (subscript 2). For 

example Ir.i = J2k,cri^ " kjvP)V'fe,<T,_R,iV'fc,fT,i?47 with 
il>k,(j,R is the operator creating a two-component Nambu 
right mover, and '4'k,a,R.,i denotes its upper component, 
rjv and (1 — \rA\) are order 1 in 1/Tk- We could work 
with these scattering states but we prefer to use the zero- 
th order scattering states, (i.e. r^v = and jr^l — 1). 
This is at the price of having to introduce an elastic 
scattering Hamiltonian. The fixed point Hamiltonian is 
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-ffj'„f where Hq is the Hamiltonian for 



a perfect scatterer (zero**^ order in 1/Tk) and 
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H' 



efc + e/c' 



i—1 k.k' ,fT 
X a2{ek + <ik')^k.ai^k','^,i 



(7) 



witliS^,,^^ - (*L,i?,» + 'i>L,L,»)/V2, where $[._,^^_, are 
the Nambu components of operators creating a scattering 
state to zero**^ order in 1/Tk- ^2 is adjusted so that H'^ 
gives the same dephasing 5^ to order 1/Tk as Hq (Eq. 
(4)). a2 is proportionnal to ai; we obtain a2(e) = ^e/A, 
with A = 2. The two-body Hamiltonian is written in the 
form 
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(8) 



[32 depends on energy e. For the same reason 
that (3i and ai were not independent, [32 and a2 
are also linked. We find /32(e) = Be/ IS. with 
B — 2. For calculating the inelastic part, it suf- 
fices to expand to second order in (32 the quantity 

(1/2) E,=±(7^c/(i") e^p(^-tJcHU{t^^)dt^^)), where rj 
and 771 are Keldysh indices, C denotes the Keldysh con- 
tour and Tc the corresponding time ordering. Lumping 
these with the elastic backscattering contribution results 
in the following expression for the averaged current 
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FIG. 3: Conductance in units of ie'^/h vs. eV/A for A/Tk = 
0.2, taking into account the many-body terms (solid line), in 
the resonant level approximation (dash-dotted line) and for 
infinite Tk (BTK result), (dashed fine). 
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As usual, this can be retrieved in a simple way by using 
Fermi golden rule for the elastic process and the three 
inelastic processes, {/R-,TVj {TZ,C), {TZ,TZ) i^,^) 
and {TZ,C) (C,^), where 7^ is a right-mover and £ 
a left-mover. For example, for the second process, two 



4 



right-movers are backscattered; one right-mover carries 
a charge e for the incoming electron and also e from the 
Andreev reflected hole, so 2e in total. The current is of 
the form (1/2) d{Nc — N-ji) / dt, where N^r ^j^-^ is the oper- 
ator number for left (right) movers and the change in the 
current is 4e. The diagram multiplicities are the same as 
for the normal case but phase space integrals are differ- 
ent because of the e dependence of a2 and (32 ■ Collecting 
the four contributions gives the average current 

(/) = 2 l\al{t){2e)T^{€)m^ 

de, (10) 



^/3|(6)e*r0,,(6) 



where the factor 2 comes from the spin; Ta{e) and 
the are 2TT/h times the right mover self-energies 

due to the above mentionned processes and m is the 
multiplicity of the diagram. For instance, for the sec- 
ond process, m2 = 2 and = 4e, so that the 
contribution is 2 i?^(7r^i'^r|.)^^ m2e* z^^ J, 



jf^yie/An2n/n) (1/4) r_^^ dx J^-^y de. 



with J ~ 
For the 

first and third processes, el = = 2e and mi = m2 = 4, 
but the self-energy triple integrals obtained by integrat- 
ing Tp_j{e) on energy e give a contribution 1/16 smaller 
than for (7^,7^) ^ (£,£). 

These results are summarized in Fig. 3, showing the 
conductance g = {I)/V versus eV/A, in units of 4e^//i, 
for A/Tk = 0.2. Many-body terms make g decrease. In 
the case oi A/Tk = 0, the BTK resulti^ is retrieved. For 
comparison, the result given by a resonant level formula- 
tion, without the many-body terms is shown. 

Noise: We now turn to the zero frequency shot-noise 
calculation. As, to zero-th order in a2 and [32, there is 
no partition noise (unlike SU(4) case^'^), applying the 



Schottky formula for each process is sufficient^. As in 
the normal case, direct inspection of the four processes 
gives 



[alie) (2e)2r„(e)m„ 



+ 5]/32'(e)(e;)2r^.,(6) 



de. 



(11) 



The resulting Fano factor with A = B = 2 is thus 10/3, 
to second order in cV/Tk for any eV < A. 

The finite frequency noise is now discussed. For a nor- 
mal NS junction, at zero temperature, the absorption 
finite frequency noise S{uj) goes to zero for lo > 2eV. 
Here, within perturbation theory in Keldysh to second 
order in this is again the case. However, inspec- 

tion of the diagrams suggests that inelastic terms bring 
more noise than in the case of an imperfect NS junction 
having the same Andreev conductance. 

In conclusion, we have proposed a generalization of 
the floating of the Kondo resonance to the case where 
one electrode is superconducting, in the regime of small 
A/Tk, close to the unitary limit. This enables a calcula- 
tion of the many-body correction to the subgap Andreev 
current and zero frequency noise. Remarkably, the Fano 
factor sticks to the value 10/3 as long as eV < A. In 
the region eV > A, not studied here, it is expected to 
decrease and to reach eventually the normal state limit 
5/3 for ey > A. 
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